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Abstract: We improve a previous quenched result for heavy- light pseudoscalar meson 
decay constants with the light quark taken to be the strange quark. A finer lattice res- 
olution (a ~ 0.05 fm) in the continuum limit extrapolation of the data computed in the 
static approximation is included. We also give further details concerning the techniques 
used in order to keep the statistical and systematic errors at large lattice sizes L/a under 
control. Our final result, obtained by combining these data with determinations of the 
decay constant for pseudoscalar mesons around the D s , follows nicely the qualitative ex- 
pectation of the 1/m-expansion with a (relative) 1/m-term of about — 0.5 GeV/mps- At 
the physical b-quark mass we obtain Fb s = 193(7) MeV, where all errors apart from the 
quenched approximation are included. 
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1. Introduction 

A big experimental progress is expected for the next years in flavour physics, mainly due 
to the new LHCb experiment [1] at CERN, CDF at Tevatron and the future super-B 
factories [2]. It will then be possible to determine with a precision of a few percent all 
entries of the CKM matrix, which describes flavour changing currents in the Standard 
Model. For this programme to be successful and to provide constraints on New Physics 
beyond the Standard Model, accurate theoretical predictions to be compared with the 
experimental results are extremely important. 

Lattice QCD is the most appropriate tool for such computations, as they involve matrix 
elements of the operators in the Weak Effective Hamiltonian among hadronic states and 
they therefore require a non-perturbative approach. 

Still, b-quarks on the lattice pose a two-scale problem: the lattice spacing a must be 
smaller than l/rrib and the size L must be large enough such that the physics is not distorted 
by finite-size effects. Heavy Quark Effective Theory (HQET) on the lattice, as formulated 
by Eichten and Hill in refs. [3, 4], allows to circumvent the problem in a theoretically sound 
way. Formally, it consists in an expansion of the QCD action and correlation functions in 
inverse powers of the heavy quark mass. 
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Numerical applications have been plagued for a long time by the exponential growth of 
the noise-to-signal ratio in static-light correlation functions. This is due to the appearance 
of power divergences in the effective theory. Such divergences are non-universal, and we 
have given in [5] first evidences how the problem can be overcome by minimally modifying 
the Eichten-Hill discretization of the static action. 

The result has been substantiated by successive studies [6, 7, 8, 9] and also in the theory 
with TVf = 2 dynamical quarks [10]. Most of the cases dealt with the non-perturbative 
renormalization of static-light operators in a finite-volume scheme. 

Here we show that also for physically large volumes fine lattice resolutions can be 
reached and precise results obtained. To emphasize the importance of such studies, we 
find that the result for Fb s in the continuum limit changes by one standard deviation (7%) 
of the result quoted in ref. [5]. At the same time, of course, we reduce the systematic 
uncertainty owing to the extrapolation to zero lattice spacing. 

The dependence of the decay constant on the heavy quark mass can very well be 
reconstructed by combining the continuum static result with continuum results at quark 
masses around the physical charm quark mass. The connection of the two different regimes 
is smooth, once the renormalization and matching of the static result is taken into account 
with sufficient precision. 

The paper is organized as follows. In Section 2 we describe the numerical setup used 
in the static approximation. Section 3 deals with the fitting procedure adopted to extract 
effective energies and matrix elements. Results in the static approximation are collected 
in Section 4, while Section 5 contains details about the simulations with relativistic quarks 
around the charm and our central results. Conclusions are drawn in Section 6. The discus- 
sion of contaminations from excited states in static-light correlation functions is deferred 
to the appendix. 

2. Heavy-light hadron physics in the static approximation 

We consider a heavy-light meson system in the framework of quenched 0(a) improved 
lattice QCD with Schrodinger functional boundary conditions, where the heavy quark 
flavour is treated at the leading order of HQET, the static approximation. The basic setup 
of the computations presented here follows our earlier determinations of hadron properties 
by means of numerical simulations of the QCD Schrodinger functional in physically large 
volumes, see e.g. refs. [11, 12, 13] for studies in the strange and charm quark sectors 
and refs. [5, 6, 7, 10] for B-physics applications. 

A particularly important technical ingredient of extracting B-meson masses and matrix 
elements from lattice HQET is the use of the alternative discretizations of the static theory, 
which were introduced in refs. [5, 6] to temper the well-known problem of exponential 
degradation of the signal-to-noise ratio encountered in static-light correlation functions 
when computed with the traditional Eichten-Hill [4] lattice action for the static quark. 

Among the static quark actions studied in detail in ref. [6], in the present work 
we restrict ourselves to the "HYP-action" , which turns out to yield the largest gain (of 
more than one order of magnitude compared to Eichten-Hill at time separations of about 
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1.5 fm) in the signal-to- noise ratios of static- light correlation functions. This action is con- 
structed from the Eichten-Hill action through replacing the temporal parallel transporters 
in the lattice derivative acting on the heavy quark field by the HYP-link that is obtained 
by a sensibly chosen smearing prescription for the gauge links located within the neigh- 
bouring hypercube [14]. Actually, in the context of static quarks, there are two favourable 
parameterizations of the HYP-link available (referred to as "HYP1" and "HYP2" in the 
following), with the second being even superior to the first. For more details the reader 
may consult ref. [6]. 1 

2.1 Correlation functions and their quantum mechanical representation 

Our starting point are B-meson correlation functions defined in the Schrodinger functional 
with a vanishing background gauge field [15]. Quarks with finite masses, also referred to as 
relativistic quarks, are labelled with an index "1", and the static ones with an index "h". 
The 0(a) improved axial vector current in static approximation is then defined as 

(Af at )o(x) = A s ^\x) + acT t 5A s ^ t (x) (2.1) 

with the (bare) unimproved current 

A s ^\x) =Ux) l0l ^ h (x) (2.2) 

and its 0(a) counterterm 

sa s ^\x) = Mxhn^ (v, + v*) Vh(^) (2.3) 

multiplied by the improvement coefficient c^ at . In order to suppress contributions from 
excited B-meson states to the correlation functions of interest, we implement wave functions 
w(x) at the boundaries of the Schrodinger functional such that an interpolating B-meson 
field is constructed in terms of the boundary quark fields (j and £ h . In this way, the 
correlation function of the static axial current, /jf* 1 *, takes the form 

fF*(x ,« H ) = -±((A?*)o(x)0( Ui )) , OH = ^^c h (y)75^(y-z)Ci(z). (2.4) 

The boundary-to-boundary correlator /f tat , 

f^\T,u i ,u j ) = -±(o'(u i )0(u j )) , o'H = |!^c;(y)75^(y-z)caz), (2.5) 

serves to cancel the overlap of the (boundary) interpolating fields O(uji) with the B-meson 
state, as made explicit in eqs. (2.10) - (2.12). 

x As the static quark action HYP2 became available only at a final stage of this project, only the 
simulations at our finest lattice resolution (/3 = 6.45) were done for both HYP-actions, HYPf and HYP2. 
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As for the choice of cc>(x) itself, we follow [5] and opt for a set of four hydrogen-like, 
spatially periodic wave functions 

Wi(x) = -^r ^ w * d x ~ nL D ' « = !,••-, 4, 
1 nez 3 

^(r) = r~ 3/2 e- r / ao , U 2 (r) = r^' 2 , 

aj 3 ( r ) = r 5/2 re- r / (2a ^ , uJ 4 (x) = L~ 3 / 2 , (2.6) 

with ao = 0.1863ro, the hadronic radius tq = 0.5 fm [16] and the normalization factors iVj 
fixed by a 3 ^ x w?(x) = 1. Apart from investigating the Schrodinger functional correlators 
for single wave functions, it is now also possible (and — as will become clear in the next 
sections — even advantageous in practice) to form suitable linear combinations of two of 
them in order to cancel the first excited state in the B-meson channel almost completely. 

For large T and xq, the correlation functions / A tat and /f tat allow for a computation 
of the pseudoscalar decay constant in the static approximation through the expression for 
the local renormalization group invariant (RGI) matrix element of the static axial current, 

c< GI (*o^) = -Z RGI (1 + bTam q ) 2L 3 / 2 e to-T,2)EMw) , (2 . 7) 



where the effective energy E e g reads 



E cS (x ,uJi) = — In 



/ A tat (x -a,^) (2g) 



2a [/^(xo + a,^) 

The 0(a) improvement coefficients c A at (cf. eq. (2.1)) and 6 A at depend on the discretization 
prescription of the static theory and have been perturbatively determined in ref. [6]. In 
eq. (2.7), Zrgi is the renormalization factor that relates the bare matrix elements of ^Q tat 
to the RGI ones. It is non-perturbatively known from ref. [17] for the relevant range of bare 
couplings employed here and has a negligible uncertainty in comparison to the statistical 
error associated with the bare matrix element. As long as <C xq <C T, the local RGI 
matrix element is expected to exhibit a plateau, from which eventually the value of ^rgi 
to enter the formula 

FpsV^ = Cps (M/A m ) x $ RGI + O (1/M) (2.9) 

for the pseudoscalar decay constant, Fpg, can be extracted. In this equation, mps is the 
meson mass, and the conversion function Cps [17, 18] translates the RGI matrix elements 
of the static effective theory to the corresponding QCD matrix elements at finite values of 
the heavy quark mass. It is parameterized in terms of the RGI mass of the heavy quark 
(M) and the QCD A-parameter in the MS scheme (A^g) [19]. 

Before coming to describe our analysis procedure for the computation of ^rgi an d the 
static binding energy based on simulation results for the correlation functions in eqs. (2.4) 
and (2.5), let us have a look at the quantum mechanical representation of these correlators. 
First, we write down the expressions for / A tat and /f tat for large values of xq and T — xq, 
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while L remains arbitrary at this stage. We neglect terms of order exp{ — (T — xo)w-g}, 
where the energy difference itlq = E[ 0) - 4 0) is the mass of the ++ glueball. In this 
approximation we obtain the decompositions [11] 

-2/l tat (*o,"i) - JZ^e-^ > & ] = 7 (fc) «f ) , (2-10) 

k>0 
fe>0 

Here, the energy Eq of the lowest state can be identified with the binding energy E stat of 
the static-light system, while in addition we have introduced 

( fc)= (fe,PS|ips(a;0) > 7 W = <0> 0|Ar|fc,PS> (2-12) 
\U,u I io / 

in terms of the k-ih excited static B-meson state, |/c,PS), the vacuum | 0, 0) and the 
boundary states [11] | ips(wj) ) and | io ), all in the finite-volume normalization (tpltp) = 1. 
The desired static-light matrix element, related to the decay constant according to eq. (2.9), 
is then encoded in 7^°) , because we have 

Zrgi (1 + bf at am q ) s/2L 3 / 2 x 7^ = Z RG1 (l + 6^ at am q ) $ bare = $ RGI . (2.13) 

From the above definitions of the correlators one infers that <J>rgi is of mass dimension 
3/2; thus, r^ 2 $RGi is dimensionless. 

2.2 Observables 

According to the foregoing discussion, the Schrodinger functional correlation functions in 
eqs. (2.4) and (2.5) obey the following asymptotic behaviour for large xq: 

-2/I tat (x ,u; l ) Xo ~°° fe-«^ + f e-^^^) , (2.14) 

2/f at (T',u W ) T '~°° a^af e~ T ' E ^ , (2.15) 

where A stat denotes the energy gap to the first excited state in the pseudoscalar channel. 
In (2.15) it is already assumed that, for our values of T', contributions to /f tat from excited 
states can be neglected — an assumption justified by a numerical analysis in appendix A. 
Above, we further exploit the freedom to choose a different time extent T' / T for the 
calculation of /f tat ; the reason of this choice will become clear in the next section. From 
the large-time asymptotics of /^ tat and ff tat one concludes that E stSu t and <1?rgi can be 
obtained from the associated asymptotic behaviour of E e g (xo, <^i) and ^^ Gl (xo,u>i) as 

„ , \ x ^oo „ sinh(aA stat ) _ a;nA stat 

Eefr xo,^ °~ ^stat + ^W -e xoA , 2.16 

# 0) a 



(0) 



$ff GI (x ,Wi) ~ $RGI < 1+ , 

I Pi 



, , sinh(aA stat ) 



(x - 



• (2-17) 
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set 


& 


a [fm] 


L/a 


T/a 


T'/a 






A, A' 


6.0 


0.093 


16 


24 


20 


0.133929 


0.135196 


B, B' 


6.1 


0.079 


24 


30 


24 


0.134439 


0.135496 


c, a 


6.2 


0.068 


24 


36 


30 


0.134832 


0.135795 


D, D' 


6.45 


0.048 


32 


48 


40 


0.135098 


0.135701 



Table 1: Simulation parameters for the calculation of the static- light correlation functions. Un- 
primed (primed) data sets refer to volumes V — L 3 x T (V = L 3 x T"), and the statistics varies 
between 0(1000) measurements for set B' and 0(2500 — 5000) measurements for the other sets. All 
simulations employed the static action HYP1, except for (3 — 6.45, where both versions, HYP1 and 
HYP2 were used. The global periodicity phase in spatial directions of the quark fields [20] is set to 
= in all cases. 

2.3 Simulation details 

The simulation parameters of our data sets are summarized in table 1. The corresponding 
quenched gauge field ensembles have been generated by a standard hybrid overrelaxation 
algorithm, where each iteration consists of one heatbath step followed by a few (in our 
case five) microcanonical reflection steps, and the sequential evaluations of the correlation 
functions were separated by at least 5 — 10 iterations. 

Thanks to the precise knowledge of the RGI strange quark mass in the quenched 
approximation, M s [12], as well as the O(a) improved relation between the renormalized 
PCAC quark mass and the subtracted bare quark mass am q = \ (k^ 1 — k^ 1 ) for the relevant 
range of bare couplings j3 = Q/qq [21], the mass of the light flavour can be directly fixed 
to the strange quark mass by proper choices of the hopping parameter, k = k s , without 
any need for interpolations in the light quark mass. More concretely, our values for k s at 
each P were obtained by solving M s = ZyiZ{\ + & m am q )m q for n, where Zm is known 
from [19], Z and b m from [21] and the critical hopping parameters from [12]. 

3. Analysis strategies for the static-light sector 

In this section we describe the extraction of the matrix element of the static-light axial 
current and of the static binding energy from the correlation functions. We start with E stat 
and display the effective energies for the data sets with (3 = 6.45 and all wave functions 
in figure 1. 

Apart from analyzing the correlation functions and the observables deriving from them 
separately for each wave function, it is advantageous to construct linear combinations, 

u v = W~ ( Ui + pi 3 ^i) ' ( 3 - x ) 

which enhance the quality and the extent of the plateau in the effective energy by (approx- 
imately) eliminating the contribution of the first excited state. Since 

Pij = -PV/PP (3-2) 
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Figure 1: Effective static-light binding energy from data set D ((3 — 6.45, 32 3 x 48) for the static 
actions HYP1 (left) and HYP2 (right). 



achieves this goal exactly, a practical way for finding such linear combinations is to perform 
a simultaneous fit to eq. (2.14) for all values of i, with /3^°\ , -Estat and A stat as fit 
parameters. Of course, the fit range has to be chosen with care. We list it together with 
the thus determined pij in table 2. Errors on the coefficients pij are shown for illustration, 
but we continue the analysis with just the central values. The so optimized correlation 
functions are subsequently analyzed without assuming anything on which excited state is 
present. In other words, we treat them just as if they were arbitrarily chosen trial wave 
functions. In the above analysis we excluded W3, since the time dependence of its effective 
energy, cf. figure 1, suggests that several states contribute in our chosen range for xq. Still, 
we checked that upon including also 103 the final results of the next section are unchanged 
within their uncertainties. 

After replacing the set of a;, with the linear combinations Uij, the static binding energy 
-E'stat is extracted by fits to eq. (2.16). Either this is done by dropping the correction term 
and fitting to a constant in a rather restricted interval, or we allow for smaller xq where 
deviations from a plateau are visible and include the correction term. These fits again may 
be performed simultaneously to all linear combinations with common fit parameters for the 
two energies. Obviously, having switched to the linear combinations of wave functions, the 
fit parameters in the correction term are expected to refer (approximately) to the second 
excited state. Indeed, for instance at j3 = 6.45, we now obtain a value of aA stat « 0.60, 
instead of aA stat ~ 0.11 for the original single wave functions with the HYP1 action. 

As for the effective energy, also for the RGI matrix element of ^Q tat we have investi- 
gated several methods of computing this quantity from the static-light correlation functions 
at hand. For all lattices, and especially at higher (3 (and correspondingly larger T/a), we 
observe the dominant part of the statistical uncertainty of this quantity to be carried by the 
correlator between the Schrodinger functional boundaries, /| tat . Hence, as already antici- 
pated in Section 2.2, we have conducted additional simulations with temporal extensions 
T < T, in order to reduce the error contribution to ^rgi from ff ta,t by calculating the 
latter on lattices with smaller time extents. Taking the crucial (3 = 6.45 data point as an 
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[ tmin/ imax/ O ] 


Pij 


6.0 


HYP1 


(1,2) 


[10, 18] 


1.0(5) 




HYP1 


(1,4) 


[10, 18] 


0.23(5) 




HYP1 


(2,4) 


[10, 18] 


-0.21(5) 


6.1 


HYP1 


(1,2) 


[13,20] 


oo 
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[13,20] 


0.21(9) 




HYP1 
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[13,20] 


-0.11(10) 
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HYP1 


(1,2) 


[13, 24] 


1.4(8) 




HYP1 
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[13, 24] 


0.23(6) 




HYP1 


(2,4) 


[13, 24] 
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-0.21(6) 




HYP2 
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[15,29] 


1.5(7) 




HYP 2 


(1,4) 


[15,29] 


0.26(5) 




HYP 2 


(2,4) 


[15,29] 


-0.17(5) 



Table 2: Coefficients of the linear combinations for all data sets. The first linear combination for 
(3 = 6.1 coincides with the second wave function. 



example, f^ at (xo,uJij) is computed on a 32 3 x T/a lattice with T = 48a and subsequently 
fitted to a two-state exponential ansatz as in eq. (2.14), i.e. 

-2/l tat (x ,^) = /^e-^'/^tat + p W e -(x -T'/2)E stat - X0 A^ } (3 3) 

whereas /f tat originates from an independent evaluation of the data set generated in a 
simulation of a volume of 32 3 x T'/a, T' = 40a. The pairs (T, T") for the remaining values 
of (5 are included in table 1. 

In the present situation, the RGI matrix element is given by eq. (2.13), where $bare is 
reconstructed as 

/? {0) 

*bareK) = — ==f ■ . (3.4) 

V fl(T, LOij,LOij) 

While the non-linear fit parameters E stat and A stat are constrained by simultaneous fits 
of /| L tat to eq. (3.3), the extracted ^bare(^jj) are not constrained. We find them all nicely 
consistent. As an alternative, we also fit ^bare( x o, ^ij) = < ^barc( < ^«j) m a restricted time 
interval and found entirely consistent values. We quote the latter as our central values. 



4. Results in the static approximation 

We follow the strategy explained in the previous section. An inspection of the plots for 
E e ff(xo,uJij) is very useful to get a first impression of the quality and the extension of the 
plateaux as well as to select reasonable fit intervals for the numerical analysis. 
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Figure 2: Top: Effective static-light binding energy from data set D (/3 = 6.45, 32 3 x 48) for 
the static actions HYP1 (left) and HYP2 (right) after construction of the linear combinations of 
wave functions, labelled W12, W14 and LO24 in the panels. The final estimates of E sigit (obtained, 
as explained before, through fits within intervals given by the vertical dotted lines) are indicated 
by the horizontal lines, where the dashed lines show the error band. Bottom: Effective binding 
energies from data sets C (left) and B (right). 



Fits are performed in the range i m i n < xq < i maX) where £ max is suggested by the 
xo-dependence of E e s(xo, u>ij) and, as already experienced in ref. [5], by the observation 
that an increase of t max beyond a certain threshold (close to 3ro) is not convenient, because 
then the statistical uncertainties become too large. Sensible values for £ m ; n may be inferred 
from the plots as well; we have kept them such that i m i n > r$ = 0.5 fm [22] holds, while the 
stability of the fit parameters under shifts of t m \ n — > t m i n — tq/2 has always been checked. 
In particular, when dropping the excited state contribution, one has to take care of the 
fitting intervals to stay safely inside the plateau region of E e f[(xo, oJij). In these cases we 
usually had t m i n « 2tq. Our results are collected in table 3. The final error for E sta _ t is 
the quadratic sum of the statistical error and the difference between the values for -Egtat 
obtained by reducing t m i n by ro/2 with t max being fixed. 

The extraction of the static decay constant from the two different fits discussed in 
Section 3 yields well compatible results. We list the ones from fits to a constant in table 3 
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Figure 3: Top: Local matrix element of the static axial current from data set D (j3 — 6.45, 32 3 x 48) 
for the static actions HYP1 (left) and HYP2 (right) after construction of the linear combinations of 
wave functions, uj\2, lou and cj24- The horizontal lines reflect the result for a 3 / 2 $bare and its error 
band (obtained through fits within intervals given by the vertical dotted lines). Bottom: Local 
matrix elements from data sets C (left) and B (right). 



p 




a 3/2 $bare(Wy) 




6.0 


0.4363(13) 


0.1976(19) 


(2,4) 


6.1 


0.3986(12) 


0.1637(27) 


(2,4) 


6.2 


0.3576(13) 


0.1295(29) 


(1,2) 


6.45 


0.2852(11) 


0.0722(18) 


(1,2) 


6.45* 


0.2564(9) 


0.0691(11) 


(1,2) 



Table 3: Results for the binding energy and the bare static-light decay constant. The entries in 
the last row refer to the HYP2 action. 



and show them for (5 = 6.1— 6.45 in figure 3. We checked the dependence of $rgi upon the 
improvement coefficient c^ at by setting the latter to its tree-level value and repeating the 
whole analysis. The outcome of this exercise is again consistent with the numbers obtained 
before. One can thus conclude that the uncertainty in c^ at does not affect our results and 
that they are expected to have all linear a-effects removed. For (3 < 6.2 and the HYP1 
action there is full agreement with the results already published in ref. [5]. 
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(a/r )* 

Figure 4: Three- and two-point continuum limit extrapolations of the RGI matrix element of the 
static axial current, represented by the red and the green line, respectively. The HYP2 result at 
a w 0.05 fm is not included in the extrapolation but only added for comparison. The continuum 
limit from the two-point fit and the HYP2 point are slightly moved to the left in the figure. 

With the new data point at (3 = 6.45 a monotonic a-dependence is observed for 
(a/ro) 2 < 0.026. Therefore, after attaching the necessary renormalization and improvement 
factors (cf. eqs. (2.13) and (3.4)), we perform the continuum extrapolation of the RGI 
matrix element of the static-light axial current (from the HYP1 data) within this range, 
see figure 4, and within a restricted range (a/ro) 2 < 0.019, dropping the data point at 
[3 = 6.1 as well. In the continuum one 

finds Ttf SaGi = 1-624(67) andr^ $ RG i = 1.577(93) 
for the three- and two-point fit, respectively. As our final estimate we quote the result of 
the three-point extrapolation to ensure stability w.r.t. statistical fluctuations and supply 
it with the statistical error of the two-point fit as the total uncertainty to try to cover a 
systematic error due to possible higher-order terms in the a-expansion. This yields 

r 3 /2 $rgi = 1.624(93) , (4.1) 

which agrees within about one standard deviation with our results for both static actions 
at the smallest lattice spacing. Given the sizeable cutoff effects, it remains desirable to 
have results at an even finer lattice resolution to gain further confidence in the continuum 
extrapolation. 

5. Decay constant at finite heavy quark mass 

This section details how we compute the pseudoscalar decay constant in the continuum 
limit of large-volume quenched QCD with Schrodinger functional boundary conditions for 
pseudoscalar meson masses in the range mps = (1.5—2.6) GeV, composed of non-degenerate 
relativistic quarks. The heavy quark mass values thus cover a significant range around the 
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physical charm region. Since it will turn out that the decay constant connects smoothly 
to the static result, ^rgi, the entire mass region mps > 1.5 GeV is covered and as an 
application, Fb s can be extracted. 

5.1 Relativistic correlation functions and observables 

We compute the decay constant, defined through the QCD matrix element between a zero- 
momentum heavy-light pseudoscalar state and the vacuum, from the time component of 
the 0(a) improved axial vector current on the lattice, 

(A I ) li (x)=A li (x) + ac A d li P(x). (5.1) 

Here, A^(x) has the form given in eq. (2.2), Section 2.1, with the static quark replaced by a 
heavy relativistic quark, and the improvement coefficient ca in front of the symmetric lattice 
derivative acting on the pseudoscalar density P{x) was non-perturbatively computed 
in [23]. 

0(a) improved correlation functions /a and f\ are defined just like / A tat and /f tat , 
but only with the standard Schrbdinger functional boundary sources (to = L04 = constant 
in eq. (2.6)). Details are as in ref. [11]. For our choice of parameters, the effective mass 
derived from /a exhibits a clear plateau already for these standard Schrbdinger functional 
boundary sources so that in this part of our calculation we can pass on adjusting wave 
functions to improve the overlap with the ground state. Hence, we readily write down the 
expressions for the effective heavy-light pseudoscalar meson mass, 

mp S (x ) = — In 

and the corresponding effective pseudoscalar decay constant 

Fps(xo) = -Z A (l + + amj) -— L= e (*o-T/2) mps (5 3) 

For large enough T and xq they equal the pseudoscalar mass and decay constant. Also the 
renormalization constant Z A is non-perturbatively known [24], and am q< i and am qjS are 
the bare subtracted valence quark masses of the heavy and the strange quark, respectively. 
The coefficient 6a, non-perturbatively tuned in ref. [25], completes the 0(a) improvement 
of our observables. 

5.2 Simulation details 

We have generated quenched gauge field ensembles for five different lattice spacings. The 
four coarser lattices have the same /3-values as the ones summarized in table 1, Section 2.3, 
and slightly different geometries (L/af x T/a = 16 x 32, 24 x 40, 24 x 48 and 32 x 64, 
respectively. In addition we have generated an ensemble with (L/a) 3 x T/a = 48 3 x 96 
and (3 = 6.7859 which, using tq = 0.5 fm [16], corresponds to a lattice spacing of a = 
0.031 fm [26]. As in the case of the simulations for the static-light observables, we employed 
a standard hybrid overrelaxation algorithm with 8 to 24 micro canonical reflection sweeps 



/a (so - a) 
_ Ia(xo + a) 



(5.2) 
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p 


6.0 


6.1 


6.2 


6.45 


6.7859 


^•meas 


non 

00U 


zUl 


zol 






Ki 


0.134108 


0.134548 


0.134959 


0.135124 


0.134739 


K2 


0.128790 


0.130750 


0.131510 


0.132690 


0.132440 


^3 


0.123010 


0.125870 


0.127470 


0.130030 


0.130253 


K4 


0.119053 


0.122490 


0.124637 


0.128131 


0.128439 


K 5 


0.115440 


0.119370 


0.122000 


0.126330 


0.126774 


K 6 


0.112320 


0.116640 


0.119680 


0.124730 


0.123571 


K 7 


0.109270 


0.113960 


1 1 7370 


0.123120 


1 1 7625 



Table 4: Summary of simulation parameters for the calculation of heavy-light correlation functions 
with relativistic quarks. K\ = n s corresponds to the bare subtracted valence quark mass of the 
strange quark, am q . s , while Ki, i — 2, . . . , 7, refer to our choices for the bare subtracted quark mass 
of the heavy flavour within the charm region. 

plus one heatbath sweep forming one iteration. Subsequent evaluations of the correlation 
functions were separated by 100 iterations for (3 = 6.0, 6.1, 6.2 and 6.45 and by 50 iterations 
for P = 6.7859. As before, we used the non-perturbatively improved Wilson quark action 
with c sw taken from [23]. 

In contrast to the hopping parameters in the static-light simulations described in Sec- 
tion 2.3, we here have determined k s at each (3 by fixing the RGI strange quark mass to 
its value found at finite lattice spacing in [12]. 2 Therefore, ki = k s in table 4 differs from 
the choice in table 1 by 0(a 3 ). The hopping parameter for the physical charm quark is 
known for the four coarser lattices [13], and through an extrapolation it was estimated at 
(3 = 6.7859. We then have guessed further hopping parameters in the vicinity of the charm 
quark value such as to yield a homogeneous covering of the region mps = (1.5 — 2.6) GeV 
with simulation points. 

5.3 Data analysis and results 

Due to uncertainties in the determination of the simulation parameters, we do not obtain 
Fps^/mps for each lattice spacing at exactly the same values of romps- Therefore, we 
first interpolate it (at fixed lattice spacing) linearly in l/(romps) to the common points 
romps = 3.768, 4.327, 4.955, 5.653, 6.211 and 6.560, which are all close to the actual 
simulation points. In a second step we estimate the continuum limit of the decay constant, 

3 /2 

r Fps^/mps, at fixed r mp S . 

Since we expect 0(a 2 ) scaling to break down for too large values of the heavy quark 
masses [28], we follow ref. [29] and exclude data points with aM > 0.64 from the dis- 
cussion of the continuum limit. Representative examples for the a-dependence are shown 
in figure 5. Since the slope of the data at small lattice spacings is not well determined, 

2 In order to determine the hopping parameter of the strange quark for the lattice with the finest reso- 
lution, we were required to extrapolate the values for the strange quark mass in [12] and the quark mass 
renormalization constant to /3 = 6.7859. Details can be found in [27]. 
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Figure 5: Lattice spacing dependence of the relativistic decay constants for two of the larger quark 
masses. 

but on the other hand we have results very close to the continuum limit itself, we take as 
our central value for the continuum limit the (weighted) average of the data at the two 
smallest lattice spacings. To account for a possible systematic error, we then added the 
difference between the average and the datum at (a/ro) 2 = 0.018 linearly to the statistical 
uncertainty of this fit to a constant. 

The outcome of the continuum extrapolations of the relativistic decay constant for the 
available six values of the quark mass within the charm region is illustrated in figure 6. It 
nicely reflects that our continuum value for the RGI matrix element of the static-light axial 
current, given in eq. (4.1), can be perfectly combined with the associated QCD estimates 
by means of a linear interpolation down to even rather low values of the heavy-light meson 
mass of about 1.5 GeV. 

Indeed, as an effective description of the mass dependence of the decay constant, we 
fit the static result and the relativistic data points to the form suggested by the HQET 
expansion, 



where we only include those points among the relativistic data into the fit that obey 
f^PS ^ ^-Ds (cf. the filled blue circles in figure 6). Here, Cps is the conversion function, 
which relates the static effective theory and QCD and has already appeared in (2.9). It 



estimated to be smaller than our statistical errors (see figure 2 in [18]). The resulting slope 
B = —1.1(2) translates into 




(5.4) 



is taken from perturbation theory in the form of [18]; its uncertainty 




B/r 



0.45(9) GeV for r = 0.5 fm . 



(5.5) 
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l 1 1 1 1 1 1 1 1 1 1 1 1 r 




0.1 0.2 0.3 

l/(r m ps ) 

Figure 6: Interpolation of the decay constant between the result in the static limit (open red star) 
and the results in relativistic QCD (blue circles). The black squares do not include the matching 
coefficient Cps- Our final quenched result at the physical point of the B s -meson, l/(rorriB s ), is 
represented by the filled red circle: r^ 2 FB sy /mB s /Cps(Mb/ Aj^g) = 1.487(55). 



We emphasize that eq. (5.4) and thus our value for the slope has to be considered an 
effective, phenomenological description and, consequently, B is an effective slope. A true 
HQET expansion can not be denned without a non-perturbative matching of the effective 
theory and QCD, i.e. a non-perturbative definition and estimate of Cps- The reason is that 
at asymptotically large quark mass, any unknown perturbative correction in Cps dominates 
over the non-perturbative 1/m-term. We refer to ref. [30] for a more thorough explanation. 

Converting the interpolation result Tq^ 2 Fb s ^/m.B s / Cps (M^ / Aj^g-) = 1.487(55) to phys- 
ical units (using roA^jg from [19] and r Mb from [7] in the evaluation of Cps as well as the 
experimental mB s -value) yields for the quenched decay constant of the B s -meson: 

F Bs = 193(7) MeV . (5.6) 

Being obtained from our effective description, this result is, however, not affected by our 
previous cautionary remarks on the HQET expansion. All that is needed for its determi- 
nation is a safe interpolation, which eq. (5.4) does represent. 

Finally, note that in an estimate of the 1/mps-correction the mass dependence in Cps, 
i.e. the anomalous dimension of the current in the effective theory, plays a numerically 
important role. This is shown by the squares in figure 6 where, as an illustration, Cps 
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has been dropped. The difference reveals that Cps accounts for about 50% of the mass 
dependence of i*ps\Aftps in the considered region. 

For other calculations of Fb b we refer to [31, 32, 33, 34, 35, 36, 37] and references 
therein. 

6. Conclusions 

We have presented a computation of the heavy-light pseudoscalar meson decay constant 
in quenched lattice QCD reaching a precision in the continuum limit of around 4%. The 
computation is founded on the 0(a) improvement and non-perturbative renormalization 
of the relativistic theory and the static approximation of HQET carried out earlier by the 
ALPHA Collaboration. Here, we have added large-volume computations dealing with the 
heavy quark both in the static approximation and in relativistic QCD with masses around 
and somewhat heavier than the charm quark's mass. 

The b-region is reached through a well controlled interpolation linear in the inverse 
of the meson mass. The final result for Fb s is nicely consistent with refs. [38, 39], where 
different strategies are applied but the same inputs are used to fix the quenched theory. 3 
The effective linear slope leads to a 0(l/mpg) correction of about 10% at the physical b- 
quark mass. The effective linear pattern is preserved for masses in the charm region, and no 
evidence is found for 0(l/mp s ) corrections at the precision level of a few percents. These 
conclusions on the mass dependence of the decay constant depend on a precise enough 
knowledge of the conversion function Cps [18] made possible through the perturbative 
result of ref. [40]. It will be very interesting to compare the present result with a direct 
HQET computation including 1/m-corrections [41]. 

At the more technical level, an important role towards a sensitive reduction of the 
statistical uncertainty is played by the choice of two different physical time extents for the 
axial current correlator and for the boundary-to-boundary correlator in the effective theory, 
as well as in the construction of interpolating fields with wave functions. By exploiting 
the latter and the fitting methods used for the decay constant we also computed the 
binding energy of the static-light system for the HYP1 action in the range of couplings 
6.0 < (3 < 6.45 and at (3 = 6.45 for the HYP2 action. They are of interest in a computation 
of the b-quark mass [42, 7, 38, 39]. 
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A. Ground state dominance for ff 



stat 



In the description of the analysis method to extract 3>r.gi in Section 3 we have implicitly 
assumed that eq. (2.15) holds and that contributions from the first and possibly higher 
excited states to /f tat in eq. (2.11) can be neglected. 

In order to arrive at a quantitative criterion in how far this assumption is justified 
for our data, let us suppress the w-dependence of the correlation functions to lighten the 
notation and write down the quantum mechanical representation of /| tat including the first 
excited state correction, taking over the notation introduced at the end of Section 2.1: 



2 ystat 



T'E S 



1 + 



a 



(i) 



o 



(0) 



-T' / A stat ) 



(A.l) 



By virtue of eqs. (2.10) - (2.12), the coefficients are related to the (3^ appearing in 
the corresponding decomposition of / A tat through 



(A.2) 



where in the second step < 0, | Aff at | fc.PS) oc Fpf' {k) ^W^ (with Fp S at ' (0) = F$f) has 
been used. Since we expect 
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(A.3) 



we get [ck^/c^ )] 2 ~ 2 ; which leads to 4 the following correction term to ^rgi 

owing to a second, higher state possibly present in /| tat : 



A/i 



/3(°) 



-T'A S 



(A.4) 



Here, (3^ and A stat are accessible through the two-state fits of / A tat discussed in the 
main text. 

After building linear combinations of the original wave functions, we find e -T'A stat _ 
O(10 -6 ) so that the magnitude of A/i is essentially driven by the ratio of the linear fit 
parameters (3^ and [3^. For all lattices and linear combinations we found A/i to be 
orders of magnitude smaller than the statistical uncertainty associated with ^rgi itself, 
thereby supporting the validity of the one-state dominance for y| tat . As already pointed 
out in Section 3, thanks to the construction of linear combinations we actually expect 
the estimated correction term to be predominantly governed by the second excited state 
instead of the first one that is written in the formulae of this appendix. 



4 Note that /f at enters $ RG i with a power of -1/2, cf. eq. (2.7). 
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